function utility = Utility(c,y,w,pr)

if pr.gam == 1,
    utility = log(c)-pr.Om./(1+pr.sig).*(y./w).^(1+pr.sig);     
else
    utility = c.^(1-pr.gam)./(1-pr.gam)-pr.Om./(1+pr.sig).*(y./w).^(1+pr.sig);
end